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We review a method for numerical simulations of lattice gauge theories at non-zero 
baryonic chemical potential we recently proposed. We first report on a test of the method 
using a solvable model and then present results for the phase structure of four flavour 
QCD. For the first time the region of chemical potential up to 1.4 Tq is explored, finding 
a first order transition line. 



1. Introduction 

The existence and nature of deconfined phase(s) of QCD at high baryon density 
is one of the open key issues of the theory of Strong Interaction Progresses in 
this field are heavily delayed by the absence of a direct method of numerical simu- 
lation based on an importance sampling Monte Carlo scheme. In order to, at least 
partially, overcome this empasse, analytical approaches in the strong coupling limit 
using the Lagrangian El or Handltoni^H^approaches, as well as other indirect 
methods, have been proposed in the last years, which in principle allow to extract 
physical informations of the behaviour of the vacuum of QCD in the portion of the 
space parameter that is not accessible to direct numerical simulations. Between the 
previously denoted indirect methods it is worthwhile to cite the Imaginary Chemical 
Potential approach and the Doubly Reweighting method lalalilalsLLiil 

Following this line we have recently proposed a different approach to the prob- 
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lem, which we hope can greatly help to improve the knowlegde in the field HI. In 
this proceeding we present a concise description of the method, some results of the 
test carried out in an exactly solvable model (3-dimensional Gross-Neveu at non 
zero baryon density in the large N limit) and finally discuss the new results we 
obtained concerning the phase structure of Q CD with four flavour of light quarks 
at finite temperature and baryon density^! 

2. The Method 

Starting from the usual discretization of QCD with staggered fermions with the 
chemical potential term introduced a Id Hasenfratz and Karsch 



S = S PG + - ^ 51 V'n^iW (Unjlpn+i ~ E^Li^n-i) 
n i—1 

+ Tj^'tpnVoin) ~ e _, ^_o,oV>ra-o) + m^-ipnlpn (1) 

n n 

we generalize the action introducing two independent parameters x and y: the 
coefficient of the antihermitian and hcrmitian terms in the temporal part of the 
fermionic action. The usual action is recovered when x — cosh(/ia) and y = sinh(/ia) 

— - 1 3 

S = S PG + m^lpnlpn + ^ 5Z {U n ,i^n+i ~ t/^^Vn-t) 

n n i—1 

+ ^ x ^2^nVo(n) (U n fitp n+ o - Ul_ Q0 1p n -(^ 
n 

+ \vYl ^nVo(n) (U nfi 1pn+0 + ^I-O.O^n-o) (2) 
n 

As in the case of Q this action suffers, in general, from the sign problem; only 
when x is real and y — iy is imaginary the model defined through (J2J can be 
simulated with standard Monte Carlo methods. Note that the case of imaginary 
chemical potential is a particular line in this subspace, defined by x = cos(fia) and 
y — sin(/ia). Studying the phase structure in the parameter space defined from (3 , 
x and y and performing an analytical continuation to real values of y we can recover 
the phase structure of the original model. 

In Ref. ^lthe reader can find an accurate discussion of what we expect for the 
phase structure in the (x, y) and (x, y) planes and the motivations for such guess; it 
is of little interest, in our opinion, to repeat that discussion here, because we want to 
concentrate on the results concerning the critical line in the (/i,T) plane. Assuming 
this attitude, let us sketch in some detail the procecdure used to extract the critical 
line in the cases under examination in the rest of the paper, i.e. Gross-Neveu model 
in three dimension and four flavours QCD. 



I 



I 
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Not to work in a three parameter space we can look at the critical line in two 
dimensional sections; in particular we can choose to work at fixed /?, in the (x,y) 
plane or at fixed x, in the (0,y) plane. 

In the first case we work at (3 smaller than the critical value of the zero density 
case and look for a critical line that starts from the y — axis moving towards 
larger values of x and y. This can be acomplished performing simulations at fixed y 
for several values of x starting from y — and repeating the procedure for several 
times at different values of y. Having determined the critical line in the (a;, y) plane 
we can fit it with an even function of y and analytically continue the line in the 
(x, y) plane. At this point the intersection of this line with the x 2 — y 2 = 1 line, that 
corresponds to real values of the chemical potential, determines the critical value 
of fia for the given j3. Repeating the whole procedure for several values of (3 will 
determine several points on the critical line in the /ia) plane; values in physical 
units can be obtained using the perturbative two-loop beta function. 

The second alternative is simply described by exchanging the role of and x. 
We work at fixed x = cosh(/ia) and look for lines in the y) plane, analitycally 
continue to the (0, y) plane and intersect them with y = sinh(/Lta) to get the critical 
values in the (0, /ia) plane, as before. 

3. Results in the Gross-Neveu model 

The Gross-Neveu model (i.e. four-Fermi model in 2+1 dimensions) is an ideal test 
bench for our approach: in the chiral limit it posseses a phase structure in the 
(/i,T) plane that is similar to the one expected for four flavour QCD and it can be 
analytically solved in the large N limit; the exact solution allows us to determine 
the critical line in the (x, y) plane for a fixed 0. In Fig. 1 the whole procedure is 
illustrated (take in mind that the variable in horizontal axis is y for the upper curve 
and y for the two lower curves). 

Repeating this procedure for others values of j3 produces the critical structure 
of the model in the (/x, (3) plane depicted in Fig. 2. In this figure the symbols stand 
for the critical points obtained with our method, the continuous line represents 
the exact results for N t = 4 and the dashed line is what can be obtained in the 
same situation performing an analytical extension starting from results obtained 
working in the Imaginary Chemical Potential approach (in all cases the analytical 
extrapolations rely on quadratic fits). It is evident a cle ar advantage in using our 
new method for increasing chemical potential values El 

4. Four Flavour QCD 

While in the Gross-Neveu model we start from data (the critical points at imaginary 
values of y) obtained solving seminumerically the gap equation, for QCD we have to 
rely uniquely on Monte Carlo simulation. For this purpose we choose to work with 
a Hybrid Monte Carlo code for staggered fermions. The lattices we have used, due 
to the fact that our method is quite expensive from a computational point of view 
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Fig. 1. Analytical continuation of the critical line from imaginary to real y and determination of 
the physical critical fi; Gross-Neveu model at /3 = 0.7 and Nt = 4. 
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Fig. 2. Critical structure of Gross-Neveu model at Nt = 4 in the /i, fi plane. 

(for any point in the (/^,T) plane we have to run the Monte Carlo to span a two 
dimensional space), are limited to 8 3 x 4 volumes at most. We choose a quark mass 
of ma — 0.05 in order to compare with existing resul ts o btained with Imagianry 
Chemical Poterntial and Doubly Reweighting methods l^l. 

For each value of x, y and fi we run around 40000 molecular dynamics trajectories 
of unit time, measuring plaquette, chiral condensate and Polyakov loop for each 
trajectory. The simulations have been performed on Linux Clusters for a total of 
100000 CPU hours. 
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We have worked, here, in the fixed x scheme, choosing 5 values of x corresponding 
to /ia of 0.2, 0.25, 0.3, 0.4 and 0.5. For each x value we perform a search of the 
critical point in f3 for a number of y values varying from 6 to 9, using a Ferrenberg- 
Swendsen procedure to precisely determine the critical values. In all cases the signal 
from the susceptibilities of the three observables we used gave identical results for 
the extimation of the critical points. The transition line was found to be a first 
order one, with a latent heat essentially independent of the y value. 

For the following stage we fitted the critical line in the (y, (3) plane, for each x, 
with a second order even polynomial: 

P(y) = Po(x) + b 2 (x)f (3) 

It is straightforward to continue the critical line to real values of y and intersect 
with the y = sinh(/xa) vertical line to obtain the corresponding critical coupling 
/3 C (/Lta). As the latent heat was constant along the line in the (y, [3) plane, we assume 
his analytical extension shares the same behaviour; hence we can affirm that the 
transitions at different /ia are still first order, as the one at /i = 0. This procedure 
is sketched in Fig. 3 for a representative value of x (as before the horizontal axis is 
y for the upper curve and y for the lower one). 
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Fig. 3. Phase diagram at fixed x in the (y, 0) plane, its continuation in the (y, fi) plane and 
determination of the physical critical point. 

Repeating this procedure for the five values of x (or fia) we have chosen we 
obtain the critical line in the (/ia, 0) plane H^l At this point we can rewrite our 
results in physical units, as reported in Fig. 4, setting the scale with the critical 
temperature at /i = and using the two loop beta function. In this figure there 
are also shown results obtained by D'Elia and Lombardo ^ using the Imaginary 
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Chemical Potential approach as well as a fit of out results with a power law: 

^ = [1-c^/TcfY (4) 

Note that, if we take seriously this ansatz for the critical line in the (/z/T c , T/T c ) 
plane, this leads to a prediction of the critical chemical potential at zero temperature 
of 1.5T C ~ 250 MeV. 

T/T c ^ — . — i — . — | — . — i — . — i — | — i — . — i — . — n 
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Fig. 4. Phase diagram in physical units. 
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